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The Palatini f(R) gravity, is able to probably explain the late time cosmic acceleration without 
the need for dark energy, is studied. In this paper, we investigate a number of f(R) gravity theories 
in Palatini formalism by means of statefinder diagnosis. We consider two types of f(R) theories: (i) 
f(R) = R+aR m -/3R~ n and (ii) f(R) = R+alnR+f3. We find that the evolutionary trajectories in 
' the s — r and q — r planes for various types of the Palatini f(R) theories reveal different evolutionary 

| properties of the universe. Additionally, we use the observational H (z) data to constrain models of 

. f{R) gravity. 
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(N ; I. INTRODUCTION 

Recently, the discovery of the acceleration of cosmological expansion at present epoch has been the most principal 
achievement of observational cosmology. Numerous cosmological observations, such as Type la Supernovae (SNIa) 
Cosmic Microwave Background (CMB) [2j and Large Scale Structure Q, strongly suggest that the universe is spatially 
i-C ■ flat with about 4% ordinary baryonic matter, 20% dark matter and 76% dark energy. The accelerated expansion of 
the present universe is attributed to the dominant component of the universe, dark energy, which not only has a large 
Q . negative pressure, but also does not cluster as ordinary matter does. In fact, it has not been detected directly and 
there is no justification for assuming that dark energy resembles known forms of matter or energy. A large body of 
recent work has focussed on understanding the nature of dark energy. However, the physical origin of dark energy as 
i well as its nature remain enigmatic at present. 

The simplest model of dark energy is the cosmological constant A [4|, whose energy density remains constant with 
y—i • time pa = A/87rG (natural units c = h = 1 is used throughout the paper) and whose equation of state (defined 
as the ratio of pressure to energy density) remains w = — 1 as the universe evolves. Unfortunately, the model is 
burdened with the well known cosmological constant problems, namely the fine-tuning problem: why is the energy of 
the vacuum so much smaller than we estimate it should be? and the cosmic coincidence problem: why is the dark 
, energy density approximately equal to the matter density today? These problems has led many researchers to try 
a different approach to the dark energy issue. Furthermore, the recent analysis of the SNIa data indicates that the 
time dependent dark energy gives better fit than the cosmological constant. Instead of assuming the equation of state 
w is a constant, some authors investigate the dynamical scenario of dark energy which is usually described by the 
dynamics of a scalar field. The most popular model among them is dubbed quintessence Q , which invoke an evolving 
scalar field <j> with a self-interaction potential V{<f>) minimally coupled to gravity. Besides, other scalar-field dark 
energy models have been studied, including phantom Q, tachyon Q, quintom ghost condensates @, etc. Also, 
there are other candidates, for example, Chaplygin gas which attempt to unify dark energy and dark matter [Toj . 
braneworld model which explain the acceleration through the fact that the general relativity is formulated in five 
dimensions instead of the usual four 

On the other hand, recently, more and more researchers have made a great deal of effort to consider modifying 
Einstein's general relativity (GR) in order to interpret an accelerated expansion of the universe avoiding the existence 
of dark energy. As is well known, there are numerous ways to generalize Einstein's theory, in which the most famous 
alternative to GR is scalar-tensor theory There are still various proposals, for example, Dvali-Gabadadze-Porrati 
(DGP) gravity [l3j . f(R) gravity [3], and so forth. The so-called f(R) gravity is a straightforward generalization 
of the Einstein-Hilbert action by including nonlinear terms in the scalar curvature. It has been shown that some of 
these additional terms can give accelerating expansion without dark energy [l5| . 

Generally, in deriving the Einstein field equations there are actually two different variational principles that one 
can apply to the Einstein-Hilbert action, namely, the metric and the Palatini approach. The choice of the variational 
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principle is usually referred to as a formalism, so one can use the metric formalism and the Palatini formalism. In 
the metric formalism, the connection is assumed to be the Christoffel symbol defined in terms of the metric and the 
action is only varied with respect to the metric. While in the latter the metric and the connection are treated as 
independent variables and one varies the action with respect to both of them. In fact, for an action which is linear 
in R, both approaches are equivalent, and the theory reduces to GR. However when the action includes nonlinear 
functions of of the Ricci scalar R, the two methods give different field equations. 

It was pointed out by Dolgov and Kawasaki that the fourth order equations in the metric formalism suffer serious 
instability problem [lfjj, however, the Palatini formalism provides second order field equations, which are free from 



the instability problem mentioned above (17|. Additionally, for the metric approach, the models of the type f(R) 



R— f3/R n are incom pat ible with the solar system experiments [18( and have the correct Newtonian limit seemed to be 
a controversial issue [19( • Another important point is that these models can not produce a standard matter-dominated 
era followed by an accelerating expansion 20, 21]. While, for the Palatini approach the models satisfy the solar system 
tests but also have the correct Newtonian limit [221 ] . Furthermore, in Ref. [23[ it has been shown that the above type 
can produce the sequence of radiation-dominated, matter-dominated and late accelerating phases. Thus, as already 
mentioned, the Palatini approach seems appealing though some issues are still of debate, for example, the instability 
problems [22|,[24j]. Here, we will concentrate on the Palatini formalism. 

In addition, since more and more cosmological models have been proposed, the problem of discriminating different 
dark energy models is now emergent. In order to solve this problem, a sensitive and robust diagnosis for dark energy 
models is necessary. As we all know, the equation of state w could probably discriminate some basic dark energy 
models, i.e., the cosmological constant A with w = — 1, the quintessence with w > — 1, the phantom with w < — 1, 
and so on. However, for some geometrical models arising from modifications to the gravitational sector of Einstein's 
theory, the equation of state w no longer plays the role of a fundamental physical quantity and the ambit of it is 
not so clear, thus it would be very useful to propose a new diagnosis to give all classes of cosmological models an 
unambiguous discrimination. In order to achieve this aim, Sahni et al. [251 ] introduced the statefinder pair {r,s}, 
where r is generated from the scalar factor a and its higher derivatives with respect to the cosmic time t, and s is 
expressed by r and the deceleration parameter q = —aa/a 2 . Therefore, the statefinder is a "geometrical" diagnostic 
in the sense that it depends upon the scalar factor and hence upon the metric describing space time. According to 
different cosmological models, clear differences for the evolutionary trajectories in the s — r plane can be found, so 
the statefinder diagnostic may possibly be used to discriminate different cosmological models. In recent works [26|, 
the statefinder diagnostic have been successfully demonstrated that it can differentiate a series of cosmological models 
, including the cosmological constant, the quintessence, the phantom, the Chaplygin gas, the holographic dark energy 
models, the interacting dark energy models, etc. 

In this paper, we focus on the f(R) theory in Palatini formalism and consider a number of varieties of f(R) theories 
recently proposed in the literature. Moreover, we apply the statefinder diagnosis to such f(R) theories. We find 
that the models in the Palatini f(R) gravity can be distinguished from dark energy models. In addition, we use the 
observational H(z) data derived from ages of the passively evolving galaxies to make a combinational constraint. 

This paper is organized as follows: In Sec. 2, we briefly review the f(R) gravity in Palatini formalism and study 
the cosmological dynamical behavior of Palatini f(R) theories. In Sec. 3, we apply the statefinder diagnosis to various 
f{R) gravity. In Sec. 4, we obtain the parameters from the observational constraints. Finally, the conclusions and the 
discussions are presented. 



II. THE PALATINI f(R) GRAVITY AND ITS COSMOLOGY 
A. a brief overview of f(R) gravity in Palatini formalism 

We firstly review the Palatini formalism from the generalized Einstein- Hilbert action 

S = Yk j d4x ^~9f{R) + SmCff^.VO, (i) 

where k = 87rG, G is the gravitational constant, g is the determinant of the metric f(R) is the general function 
of the generalized Ricci scalar R = g^ v 'i?^„(r A and S m is the matter action which depends only on the metric g^ v 
and the matter fields ip and not on the independent connection r A ^„ differentiated from the Levi-Civita connection 
{ x t j ll y}- In our paper, nature unit c = 1 is used. It should be noted that GR will come about when f(R) = R. 
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Varying the action with respect to the metric g^ v and the connection T x ^ v respectively yields 

f{R)R^ - \f{R)g„u = kT^, (2) 
v A (V=5/' W) = o, (3) 

where f'{R) = df /dR, Va denotes the covariant derivative associated with the independent connection T x ^ u and T^ v 
is the energy-momentum tensor given by 



2 5S n 



If we consider a perfect fluid, then T Ml/ = (p + p)u^u v + pg^ v , where p and p respectively denotes the energy density 
and the pressure of the fluid, is the fluid four- velocity. Note that T — g^ v T^ u = —p + 3p. 
According to Eq. (3), we can define a metric conformed to as 

h>nv = f'(R)g^- (5) 

Then, we can get the connection T x ^ v in terms of the conformal metric h^ v 

r fiv h (han^is hfjj^ ^ (6) 

furthermore, it can be equally written as 

rV = j^g Xa [Mf'(R)9,a) + W(R)g^) - dM'(R)g^)}- (7) 

Meanwhile, the generalized Ricci tensor is 

r> pA pA _r_ pA per pA per /q\ 

and thus, we can rewrite the generalized Ricci tensor expressed by the Ricci tensor R fiu (g) associated with g^ v as 

3 1 

2 Tp) 5 

-J^ V ^^'^-lj^f^^ a f'( R )' (9) 
where V M is the covariant derivative defined with the Levi-Civita connection of the metric. 



iV = RM + oTTTW (V M /'(i?))(V,/'(i?)) 



B. FRW cosmology of the Palatini /(i£) gravity and numerical results 

Since measurements of CMB suggest that our universe is spatially flat [13, we start our work with a flat 
Friedmann-Robertson- Walker (FRW) universe with metric 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (fO) 

where a(t) is the scalar factor and t is the cosmic time. 

As a result, making use of Eqs. (2) and (9), the modified Friedmann equation can be derived as 

lf(R) a l K(p + 3p) + f(R) 
1 2f'{Ry 6 f'(R) ' { 1 

where H is the Hubble parameter and the dot denotes the differentiation with respect to the cosmic time t. 
In addition, taking the trace of Eq. (2), gives 

f'(R)R-2f(R) = K T. (12) 

If we consider the universe only containing dust-like matter, then T — ~p m , where, p m denotes the energy density of 
matter. Furthermore, combining Eq. (12) with the energy conservation equation of matter p rn + 3Hp m = 0, we can 
express R as 

• 3H(f'(R)R-2f(R)) 
f"(R)R-f'(R) ' 
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where f"(R) = d 2 f/dR 2 . Replacing R in Eq. (11) with Eq. (13), yields 



H 2 _ 1 3/(fl) - f'(R)R 



6f'(R) [i _ 3 f"(R)(f'(R)R-2f(R)) ]2 : 



2 f'(R)(f"(R)R-f'(R))- 

from which, we can get the modified Fricdmann equation with respect to R given the form of f{R). 

Using the redshift z = - — 1 (usually, the scalar factor a of today is defined do = 1), the expression p m = /9 m o(l + z) 3 
and 4| = — H(l + z), we can rewrite Eqs. (12) and (13) as 

f'{R)R - 2f(R) = -3H 2 Q m0 {l + z) 3 , (15) 
dR _ 9H$Q m0 (l + zf 



dz f"(R)R-f'(R) 



(16) 



where Q m o = K ^jp° and the subscript throughout the paper denotes the present time. Therefore, we can express the 
Hubble parameter H in terms of z as 

H 2 _ 1 3fi m o(l + zf + f(R)/H 2 



#o 6/'(i?) m , 9 ggn„o(i+s)»/"(fl) p ' 

° Jy ' i L+ 2(f"(R)R-f'(R))f'(R)\ 

if the form of f(R) with respect to R is given. In order to study the cosmological evolution by using Eqs. (15)- (17), it 
is necessary to give the initial conditions: (Rq, Ho, fl m o)- From Eqs. (15) and (17), choosing units so that Hq = 1 [29j |. 
we can solve for Rq given that the rest ones except one parameter are fixed. 

On the other hand, in order to understand the cosmological evolution behavior, it is useful to define the effective 
equation of state 

Weff = -l + \{l+z)tL. (18) 

Thus, we can show the effective equation of state as a function of redshift z for any f(R) from Eqs. (16)-(18). 
In our paper, we adopt two types of f(R) theories, recently considered in the literature. 



1. f(R) theories with power-law term 

We consider the following general form for f(R) 

f(R) = R + aR m - /3R- n , (19) 

where m and n are real constants with the same sign. Such theories have been considered with the hope of explaining 
the early and the late accelerated expansion of our universe [3(1 HH . Note that not all combinations of m and n are 
agreement with a flat universe with the early matter dominated era followed by an accelerated expansion at late times. 
At the early times of matter dominated, the universe is better described by GR in order to avoid confliction with 
early-time physics such as Big Bang Nucleonsynthesis (SSN) and CMB. This implies that the modified Lagrangian 
should recover the standard GR Lagrangian for large R, and hence we demand that m < 1 and n > — 1. Now we use 
the values of a and f3 to classify such theories. 

Case 1 a ^0,/3 = 

In this case, (19) reduces to 

f(R) = R- f3R- n . (20) 

Substituting the form of (20) into Eqs. (16)-(18), with the present fractional matter density f2 m o = 0.27, the changing 
of the Ricci curvature R and the effective equation of state w e ff with the redshift z are plotted in Fig.l. It is to 
be noted that the special case of (a, n) = (—4.38, 0) corresponds to the ACDM model. We can easily see that the 
curvature and the effective equation of state decrease with the evolution of the universe for any choice of n. Moreover, 
the smaller n, the faster R decrease, and the larger the present value of w e ff. Also, the universe turns to an accelerated 
phase from a decelerated era, and tends to a de Sitter phase in the future. 
Case 2 a, /3^0 
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FIG. 2: The figures are for the model f(R) = R + ai? m — /3R~ n with m = n = 1/2. The present fractional matter density 
fimo = 0.27. The left figure is the diagram of Ricci scalar R as a function of redshift z. The right figure is the evolution 
trajectories of w e ff. 



In this case, we consider the general form of this type. To study the cosmological behavior of such theories, we 
adopt the model 

f{R) = R + aR 1 ' 2 - PR- 1 ' 2 . (21) 

From Fig. 2, we clearly find that the curvature and the effective equation of state decrease with the evolution of the 
universe for any choice of /3. the smaller j3 results in the faster decrease of R and the larger w e ff at the present time. 
It is also obvious that the universe evolves from deceleration to acceleration, and enters to a de Sitter phase in the 
future. 




FIG. 3: The figures are for the model f(R) = R + ctlnR + .The left figure is the diagram of Ricci scalar R as a function of 
redshift z. The right figure is the evolution trajectories of W e ff. 
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2. f{R) theones with logarithm term 

For this type of f(R) theory, we adopt the form 

f{R) = R + alnR + P, (22) 

which has been studied in [jjj [H| . It has been claimed that such theories have a well-defined Newtonian limit [HJ ■ 
Note that, the asymptotic behavior lim^^ +oc f(R) —> R is obtained for any choice of a and (3, and thus, the 
arbitrary a and f3 can satisfy the assumption that the universe is described by GR at the early time. However, not 
all combinations of a and (3 can explain a late-time accelerated expansion of the universe. Therefore, for the sake 
of compatibility with the observational constraints obtained in Sec. 4, we select a series values of (3, which can well 
explain the evolvement of the universe from an early-time deceleration to a late-time acceleration (see Fig. 3). 

Substituting the form of (22) into Eqs. (16)-(18), with £l m o = 0.27, the changing of the Ricci curvature R and 
the effective equation of state w e ft with the redshift z for this model are plotted in Fig. 3. Obviously, R and w e tf 
decrease with the evolution of the universe for any choice of (3. Also, the larger (3 results in the faster decrease of 
R and the larger w e tf at the present time. It is noting that, similar to the result of the above type of theories, the 
universe evolves from deceleration to acceleration, and enters to a de Sitter phase in the future. 



III. STATEFINDER DIAGNOSIS FOR THE PALATINI f(R) GRAVITY 

In this section, we turn our attention to the statefinder diagnosis. As we know, two famous geometrical variables 
characterizing the expansion history of the universe are the Hubble parameter H = a/a describing the expansion rate 
of the universe and the deceleration parameter q = —a/aH 2 characterizing the rate of acceleration/deceleration of 
the expanding universe. It is clear that they only depend on the scalar factor a and its derivatives with respect to 
t, i.e., a and a. However, as the enhancing of cosmological models and the remarkable increase in the accuracy of 
cosmological observational data, these variables are no longer to be a perfect choice. This can be easily seen from the 
fact that many cosmological models correspond to the same current value of q. As a result, the so-called statefinder 
diagnosis was introduced in order to discriminate more and more cosmological models. 

The statefinder diagnosis is constructed from the scalar a and its derivatives up to the third order. Namely, the 
statefinder pair {r, s} is defined as 

a r — 1 , , 

(23) 



aH 3 ' 3(9-1/2)' 

Since different cosmological models exhibit qualitatively different trajectories of evolution in the s — r plane, the 
statefinder diagnosis is a good tool to distinguish cosmological models. The remarkable property is that the statefinder 
pair {r, s} — (0, 1) corresponds to the ACDM model. We can clearly identify the "distance" from a given cosmological 
model to ACDM model in the s — r plane, such as the quintessence, the phantom, the Chaplygin gas, the holographic 
dark energy models, the interacting dark energy models, and so forth, which have been shown in the literatures 26j . 
Particularly, the current values of the parameters s and r in these diagrams can provide a consider way to measure 
the "distance" from a given model to ACDM model. Generally, according to the reexpression of the deceleration 
parameter q 

q = -l + (l + z)^, (24) 

where H 1 = 4^, we can also rewrite the statefinder pair {r, s} in terms of the Hubble parameter H and its first and 
second derivatives H' and H" with respect to the redshift z as 

r = l-2(l + z)^ + (l + z) 2 ^ T +{l + z fE-, (25) 
. _ -2(1 + z)H'/H + (1 + zf{H'/Hf + (1 +zfH"/H _ (2g) 



3[(1 + z)H'/H - 3/2] 

In what follows, we will apply the statefinder diagnosis to the f(R) theories mentioned in Sec. 2. The values of 
parameters we select in such theories are as same as those in the previous section. By describing the evolution 
trajectories of statefinder parameters r and s for the Palatini f(R) theories, we can differentiate f(R) models from 
dark energy models, and even discriminate various types of Palatini f(R) theories from each other. 
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FIG. 4: r(z) and s(«) for the models f(R) = R~(3R~ n , f(R) = R+aR 1/2 -/3R- 1/2 and f(R) = R+alnR+f3 with fi m0 = 0.27. 



We plot the statefinder parameters r(z) and s(z) for the above two types of Palatini f(R) theories in Fig. 4 with 
!l m o = 0.27. Fig. 4 show that r(z) and s(z) for the three models exhibit different features. For the model f(R) = 
R — (3R~ n (see Fig. 4a), the curves cross the ACDM line (r ~ 1 or s = 0) twice times and approach it in the future. 
Moreover, the curves with n > start from the region r > 1, s < 0, while the other way round, the curves with n < 
start from the region r < 1, s > 0. For the model f(R) = R + aR 1 / 2 — fiRr 1 / 2 (see Fig. 4b), the curves start from 
the region r > 1, s < 0, cross the ACDM line and tend to it. While for the model f(R) = R + alnR + f3 (see Fig. 4c), 
the curves generally do not cross the ACDM line but approach the line in the future. When the parameters a and (3 
have the same sign, the curves lie in the region r > 1, s < 0, and in reverse the curves lie in the region r < 1, s > 
for the case that a and j3 have the opposite sign. From Fig. 4, we can clearly find that the behavior of r(z) and s(z) 
for the models in the Palatini f(R) theories is different from ACDM model and even from other dark energy models. 
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FIG. 5: The evolution trajectories of r(s) and r(q) for the models f(R) = R - (3R~ n , f(R) = R + aR 1/2 - /3R~ 1/2 and 
f(R) =R+ alnR + with fl m0 = 0.27. 



Furthermore, various types of Palatini f(R) theories can be distinguished from each other in the z — r and z — s 
panels. 

The evolutionary trajectories of the statefinder pairs {r, s} and {r, q} can also help to differentiate various models in 
the Palatini f(R) theories from dark energy models. The trajectories r(s) and r(q) for the models under consideration 
are plotted in Fig. 5. It is easy to see from the figure that various models in the palatini f(R) theories in the s — r and 
q — r planes exhibit the significant differences. Furthermore, the evolutionary trajectories r(s) and r(q) also effectively 
distinguish the f(R) models from the ACDM model and other dark energy models. 
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IV. NUMERICAL ANALYSIS FROM OBSERVATIONAL DATA 



In order to make f(R) models to be compatible with cosmological observations, we now turn to constrain the 
parameters in f(R) models with the observational H(z) data. 

The Hubble parameter H{z) data depends on the differential ages of the universe as a function of redshift z in the 
form 

fl w - (27) 

which provides a direct measurement for H(z) through a determination of dz/dt (34[. By using the differential age s 
of passively evolving galaxies determined from the Gemini Deep Deep Survey (GDDS) 35] and archival data 36], 
Simon ct al. determined H(z) in the range < z < 1.8 and used them to constrain the dark energy potential and 
its redshift dependence [34) ]. In order to impose constraints on the models of f(R) gravity, we determine the best fit 
values for the model parameters by minimizing 

^ 2 _ ^ - \- H th{z<,\s) - H obs (zi)] 2 



1=1 



where H t h(zi\s) is the theoretical Hubble parameter at redshift Zi given by (17); H a t s (zi) are the values of the 
Hubble parameter obtained from the data selected by [34j (SVJ05) and cr(zi) is the uncertainty for each of the nine 
determinations of H(z). 



A. The Type f(R) = R+ aR m - /3R~ 



1. The Case of f(R) = R - /3R~ 



2. The Case of f(R) = R + aR m - /3R~ n with m = n = 1/2 



B. The Type f(R) = R + alnR + [3 



V. CONCLUSIONS AND DISCUSSIONS 



In summary, we investigate the f(R) gravity in Palatini formalism by means of the statefmder diagnosis in this 
paper. Differences of the evolutionary trajectories in the s — r plane among a series of f(R) models have been found. 
Therefore, the statefmder parameters are powerful to discriminate the models in the Palatini f(R) gravity from dark 
energy models and even from other f(R) models. Also, the q — r plane has been widely used for discussion on the 
evolutionary property of the universe. We find that the f(R) models under consideration exhibit different properties in 
the q — r plane. Since the model parameters are found to be sensitive to the f(R) models, constraining the parameters 
in the models exactly becomes a valuable task. We use the observational H(z) data to make a combinational constraint. 



Acknowledgements 

We are grateful to Xin-juan Zhang for helpful discussions. This work was supported by the National Science 
Foundation of China (Grants No.10473002, 10533010), 2009CB24901, the Scientific Research Foundation for the 
Returned Overseas Chinese Scholars, State Education Ministry. 

Note added. — As we were completing this paper we became aware of the work reported in (37| , which uses a similar 
method to discriminate the ACDM from viable models of f(R). While the authors of that work consider the models 
which is different from ours. Furthermore, we derive the statefmder parameters directly from the Hubble parameter 
H(z). In addition, we impose constraints on the models by using different observational data. 
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